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ABSTRACT 

We have obtained Spitzer IRS 5.5-35 /im spectroscopy of the debris disk around (3 Pictoris. In 
addition to the 10 /jm silicate emission feature originally observed from the ground, we also detect 
the crystalline silicate emission bands at 28 /im and 33.5 /im. This is the first time that the silicate 
bands at wavelengths longer than 10 /im have ever been seen in the /3 Pictoris disk. The observed 
dust emission is well reproduced by a dust model consisting of fluffy cometary and crystalline olivine 
aggregates. We searched for line emission from molecular hydrogen and atomic [S I], Fe II, and Si II 
gas but detected none. We place a 3a upper limit of < 17M ffi on the H 2 S(l) gas mass, assuming 
an excitation temperature of T ox = 100 K. This suggests that there is less gas in this system than is 
required to form the envelope of Jupiter. We hypothesize that some of the atomic Na I gas observed 
in Keplerian rotation around (3 Pictoris may be produced by photon-stimulated desorption from 
circumstcllar dust grains. 



Subject headings: circumstellar matter — planetary systems: formation- 
planetary disks — stars: individual (/3 Pictoris) 



planetary systems: proto- 



1. INTRODUCTION 

The A5V star (3 Pictoris (HD 39060 = HR 2020) at 
a distance d k, 19.3 pc possesses a spectacular edge-on 
debris disk imaged in scattered light a nd thermal emis- 
sion that extends to rad i i >1400AU (ISmith fe Terriiel 
fl98l Holland et aDItlM iGolimowski et all 120061) . ~An 
age of ^12 Myr has been estimated for for the central 
star based on the pro perties of late- type mem bers of the 
(3 Pic moving group (|Zuckerman et aLll2001l ). The dust 
in this system is believed to be replenished from a reser- 
voir, such as collisions between parent bodies or sublima- 
tion of comets because the estimated lifetime for the dust 
under radiation pressure, Poynting-Robertson drag, and 
grain-grain collisions is a few orders of magnitude shorte r 
than the age of the system (Backman & Parcscc 1993). 
There may be evidence for the presence of large bod- 
ies in the (3 Pictoris disk that perturb dust grains and 
parent bodies and instigate collisional cascades between 
them. A 5° warp at radii <80 AU has been observed in 
scattered light using STIS on HST and may be gener- 
ated by either a brown dwarf close to the star (at dis- 
tances <3 AU) or a 0. 17 Mj up planet at larger distances 
from the star (50 AU) (jHeatTet al.ll2000f) . The f3 Pictoris 
disk may also possess infalling, sublimating, refractory 
bodies. Time- variable, high velocity, non-periodic, red- 
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shifted atomic absorption features have been observed 
toward (3 Pictoris at ultraviolet (UV) and visual wave- 
lengths that vary on timescales as short as hours. The 
velocity of the atoms, typically 100-400 kms -1 , is close 
to the free fall velocity at a few stellar radii, suggesting 
that the absorption is produced as stellar photons pass 
through the co mae of infalling bodies at distances <6 AU 
from the star (|Vidal-Madiar et al.lH99l . 

The origin of the micron-sized dust grains required 
to account for the observed scattered light and thermal 
emission is currently not well-constrained. The broad 
9.7 /im silicate emission feature and the narrow 11.3 /xm 
crystalline olivine emission feature observed toward (3 
Pictoris appear grossly similar to those observed to- 
ward comets Halley , Bradford 1987s, and Levy 1990 XX 
(|Knacke et al.lll993l ). suggesting that the grains may be 
cometary. Models of cometary grains, idealized as fluffy 
aggregates of interstellar silicate cores with carbonaceous 
organic refractory mantles (and additional ice mantles at 
distances larger than the snow-line), in radiative equi- 
librium with the central star are able to reproduce the 
observed 10 /im silicate feature and the spectral energy 
distribution (SEP) at I R through millimeter wavelengths 
(|Li fe Greenbergj["l998f ). Spatially-resolved studies of the 
silicate emission feature, obtained using COMICS on the 
Subaru Telescope, suggest that the shape of the 10 /im 
feature changes as a function of position in the disk, with 
large and crystalline grains concentrated at the center 
and small grains located in ring- like bands at 6.4, 16 , 
and 29 AU from the central star (jOkamoto et al.1 12001 ). 
The rings of small grains may be explained by collisions 
between large orbiting bodies. High-resolution, mid-IR 
imaging of (3 Pictoris, obtained with TReCs on the Gem- 
ini South Telescope, has resolved a prominent clump on 
the south-west side of the disk at 52 AU from the central 
star which may also indicate t hat a recent collisio n has 
occured in the j3 Pictoris disk (Te lesco et al .1120051 ). 

Spatially resolved visual spectroscopy of (3 Pictoris has 
revealed a disk of atomic gas in Keplerian rotation, ob- 
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served via scattered emission from Fe I, Na I, Ca II, 
Ni I, Ni II, Ti I, Ti II, Cr I, and Cr II. The atomic 
gas possesses a NE/SW brightness asymmetry and an 
inner warp similar to that observed in the dust disk al- 
though the gas extend s to larger heights than the dust 
(Brandekcr ct al. 2004). Estimates of the radiation pres- 
sure acting on Fe and Na atoms suggest that these species 
should be not be in Keplerian rotation but should be 
accelerat ed to terminal outflow velocities ~100s-1000 s 
km/sec ([Brandeker et all 120041 ) . lLagrange et~aTI 1)19981 ) 
showed that a ring of neutral hydrogen at distance 0.5 
AU could slow d own t he radial flow of gas. However, 
iFernandez et alJ (|2006l ) have recently proposed that the 
gas will not be radially accelerated if the carbon is over- 
abundant; their model does not require a population of 
undetected atomic hydrogen. Since carbon does not ex- 
perience strong radiation pressure (-Fi-ad/^grav = /3 C arbon 
w 0) and also has a large ionization fraction (~0.5), 
they suggest that Coulomb interactions between C II and 
other ions reduce the effective radiation pressure on the 
bulk gas. In this case, the ions and neutral atoms in the 
disk couple together into a fluid, with an effective radia- 
tion pressure coefficient, that is bound to the system and 
that brakes the gas if /3 e ff < 0.5. In particular, they sug- 
gest that if the carbon abundance is >10x solar, then all 
the atomic gas will be retained. Measurements of the col- 
umn density of the stable component of atomic carbon 
(with zero velocity re lative to the star) v ia absorption 
from C I ( 3 P) A1613 (|Roberge et all 12000( 1 and absorp- 
tion from C IIA1036 and C II* A1037 superimposed on 
chromospheric O VI A 1038, suggest that the bulk of the 
atomic gas is composed of carbon with a C/Fe abundance 
ratio that is 1 6 times the solar val ue and an ionization 
fraction of 0.5 ([Roberge et al.ll2006f l. 

We report the results of a Spitzer IRS study of the 
dust and gas around f3 Pictoris, building on the model 
for the composition and sp atial distribution of the grains 
by iLi fe Greenbergj (|1998l ) and placing 3er upper limits 
on the line emission from H 2 S(2), S(l), S(0) and [S I], 
Fe II, and Si II gas. Based upon the similarity in the 
spatial distribution of the dust and gas observed in scat- 
tered light, we hypothesize that the dust and gas in this 
system are physically associated and that the observed 
gas is secondary; it has been produced from circumstel- 
lar material since the star formed. We quantitatively 
show that the observed Na I, seen in Keplerian rotation 
around the star, can be produced via photon-stimulated 
desorption in a time that is shorter than the residence 
time of the gas in the disk. 

2. OBSERVATIONS 

We obtained IR Spectrograph (IRS) (jHouck et al.1 
2004) spectral mapping mode observations of (3 Pictoris 
using the Short-Low (SL; 5.2-14 /mi; A/AA ~ 90), 
Short-High (SH; 9.9-19.6 /mi; A/AA ~ 600), and Long- 
High (LH; 18.7-37.2 /urn ; A/AA ~ 600) mod ules on the 
Spitzer Space Telescope ([Werner et al.ll2004D . The SL2 
slit has a size 3.6" x 57"; the SL1 slit has a size 3.7" x 57"; 
the SH slit has a size 4.7"xll.3"; the LH slit has a size 
ll.l"x22.3". Both the low and high resolution observa- 
tions were made with the spectrograph long slits aligned 
along the position angle of (3 Pictoris disk to within 5°. 

The SL2 (AOR key: 8972288) and SL1 (AOR key: 
8972544) spectra were obtained on 16 November 2004 



by stepping the slit in 7 and 11 1.8" intervals perpen- 
dicular to the disk, respectively, with the center position 
centered on the disk. Each SL2 and SL1 map position 
was observed using 9 and 14 cycles, respectively, of 6 sec- 
ond ramps. The SH (AOR key: 4879616) and LH (AOR 
key: 4876800) spectra were obtained on 15 December 
2003 and 4 March 2004, respectively, by stepping the slit 
in 3 2.4" and 4.8" intervals perpendicular to the disk, re- 
spectively, with the center position centered on the disk. 
Each SH and LH map position was observed using 25 
and 24 cycles, respectively, of 30 and 14 second ramps. 
High accuracy peak-up, with the blue (13.3-18.7/xm) ar- 
ray, on the nearby star HD 38891 (located 12' north of 
(3 Pictoris), was used to center the disk in the SL2, SL1, 
and SH slits; high accuracy peak-up, with the red array 
(18.5-26 /mi), on HD 38891 was used to center the disk 
in the LH slit. 

The raw data were processed into calibrated (flat- 
fielded, stray light corrected) 2-dimensional spectra with 
version Sll.0.2 of the SSC pipeline. We removed masked 
and rogue pixels from the basic calibrated data by inter- 
polating between their nearest neighbors, coadded the 
resulting images, and extracted 1-dimen sional spectra 
using the IRS team's SMART program ([Higdon et al.1 
2004) . Since separate nod observations of the sky were 
not made when f3 Pictoris was observed, we approxi- 
mated the SL sky background as a constant at each wave- 
length and subtracted this value from each pixel in the 
spatial direction across the source. Since the SL slits are 
long enough to sample the PSF in the spatial direction, 
we fit 1-dimensional Gaussians to the SL spectra in the 
spatial direction to determine the positions of the source 
in the slits. We extracted our low resolution spectra by 
summing over a window whose width varied as a func- 
tion of wavelength to account for the changing size of the 
point spread function and was centered on the source po- 
sition. Finally, we multiplied the resulting spectra with 
an a Lac relative spectral response function (RSRF), ob- 
tained by dividing an a Lac te mplate spectrum with an 
observed a Lac SL spectrum ([Jura etall[200l . Since 
point sources fill the high resolution slits, we extracted 
our high resolution spectra (without background subtrac- 
tion) by summing over all of the pixels in the SH and LH 
slits. The extracted high resolution spectra were multi- 
plied by a £ Dra high resolution RSRF. 

We constructed a spectrum that is consistent with the 
LH slit, the largest slit we used to observe j3 Pictoris, 
centered on the disk. Since the source is extended in our 
IRS observations, we added multiple SL measurements 
together to construct an aperture equal to the LH aper- 
ture. Since our SH map was not large enough to cover 
the LH central slit position, we scaled our SH spectrum 
by a factor of 1.068 to match the SL and LH spectra, 
using the LL spectrum at 5.1939 /tm to 9.9232 /mi, the 
SH spectrum at 9.9461 /tm to 19.3179 /tm, and the LH 
spectrum at 19.2636 /mi to 35.9892 /tm. The final spec- 
trum of j3 Pic (shown in Figure [T]) is that of a ll"x22" 
region centered on the star and oriented along the disk. 

The signahnoise in our spectrum is very high. Over 
spans of many spectral resolution elements, the sensitiv- 
ity of our spectrum is limited by flat-fielding errors to 
1-2% of the continuum flux density. Over smaller spans, 
< 10 spectral resolution elements or so, the sensitivity is 
limited only by photon noise from the background and 
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target, and is significantly better: 0.1-0.5% of the con- 
tinuum flux density in SH and LH, and 0.05-0.2% in SL, 
consistent with the point-to-point scatter in the spec- 
trum over such small spans, and (perforce) consistent 
with the Spitzer-IRS facility sensitivity tool, SPEC-PET. 
Thus, every spectral feature visible in our (3 Pic spec- 
trum that is at least as wide as the spectral resolution 
is real, and the upper limits on weaker spectral lines are 
significantly smaller than 1% of the continuum. For sim- 
plicity we have calculated such upper limits directly from 
the point-to-point variation of the continuum-subtracted 
spectrum. 

3. DUST PROPERTIES 

Figure [1] plots the 5.2-36 /zm IRS spectrum and the 
PHOENIX model stellar atmospheric spectrum with ef- 
fective temperature = 8000 K and surface gravity 
log g = 4.2. Figure Q] clearly shows that the dust excess 
emission above the stellar atmospheric radiation is even 
seen at wavelengths as short as A ~ 5-8 /zm, indicating 
that there may exist a population of hot dust (~ 400- 
600 K) at a distance of <10 AU fro m the central star [see 
Figure 1 of iLi fc GreenbeTd" (fl998h ]. 

Mostly prominent in Figure[T]are the broad 10 /im fea- 
ture (which reveals the presence of both amorphous sil- 
icates and crystalline silicates in the (3 Pictoris disk), 
and the 28 /jm and 33.5 /jm crystalline olivine emission 
features (s ee the inset in Figure [Ik) which closely re- 
semble the ILi fc Greenberg! (|1998l ) model spectrum (see 
their Figure 6c). 10 This is the first time that the silicate 
bands at wavelengths longer than 10 /im have ever been 
detected in the j3 Pictoris disk. 

For comparison, we also show i n Figure [1] the IRAS 
broadba nd photometry of iGillett (1986), the photom- 
etry of ITelesco fc Knackd (|1991l ) at 8.8, 10.3, 11.7 
and 12.5 /im obtained using the Big Mac spectral fil- 
ters at the NASA Infrared Telescope Facility (IRTF), 
the IRTF N band (10.1 /zm) and Q band (20 /im) 
photometry of iBackman et all (|1992l ). the IRTF 2.6- 
13.5 /im int ermediate resolution (A/AA fts 50) spec- 
trometry of iKnacke et al.l (|1993l ) . and the Kuiper Air- 
bo rne Observato r y (KA O) 47 and 95 /im photometry 
of lHarvev et alj (119981) . While the photometry of 
ITelesco fc Knackd (119911) a nd the IRTF N band pho- 
tometry of Backm an et al.l <l 19921 ) agree with the spec- 
trometry of IKnacke et al.l ([19931) very well, they are be- 
low the IRS spectrum by a factor of ~1.94 (see Fig- 
ure [TJa and Figured). This difference is much larger 
than the calibration uncertainty of IRS and is attributed 
to the fact that the Spitzer apertures (3.6" x 57" for 
SL, 4.7" x 11.3" for SH, 11.1" x 22.3" for LH) are 

10 Crystalline pyroxene Mg^Fe^^) Si03 grains also have a rich 
set of features in the mid-IR. But their peak wavelengths and rela- 
tive feature strengths are inconsistent with the IRS spectrum (e.g. 
see Chihara et al. 2001). Olivine Mg2 :E Fe2(i_ ;E )Si04 is a mixture 
of forsterite Mg2Si04 and fayalite Fe2Si04 with a mixing propor- 
tion x. In this work as well as in Li & Greenberg (1998) we adopt 
the refractive indices of Mgi.8Feo.2Si04 (Mukai & Koike 1990). 
Although in principle it is possible to infer the Mg:Fe ratio from 
the locations and strengths of the IR features (e.g. see Koike et al. 
1993, 2003, Fabian et al. 2001), neither the IRS spectrum reported 
here reveals a large number of crystalline silicate features nor does 
it cover a sufficiently broad wavelength range for constraining the 
exact crystal composition. As shown in Figure 2, the dust model 
consisting of crystalline Mgi.sFeo.2Si04 olivine reproduces the IR 
features reasonably well. 



much larger than the apertures of ground-ba s ed ob - 
servations [e.g. 3.7" diameter of Knacke et al.1 (|1993h . 
3. 9" x 4.2" of ITelesco fc Knackd ([1991|), 7.8" diameter 
of IBackman et afT i 19921 )]. and therefore Spitzer detects 
more flux from extended material than in previous ob- 
servations. 11 

The nature (e.g. size, composition, and morphology) 
of the dust and its spatial distribution in protoplan- 
etary and debris disks is mostly revealed through 
the interaction of the dust with the electromagnetic 
radiation of the central star: absorbing, scattering and 
polarizing the stellar radiation and re-radiating the 
absorbed UV/visible photons at longer wavelengths, 
ranging from the IR to submillimctcr and millimeter. 
Since the disk around (3 Pictoris subtends more than 
100", it has been studied in great detail. There exists 
a vast variet y of observa t ional data for the dust in 
this disk (see iMann et al.l (|2006f )). incl uding (1) imag- 
ing observations of scattered light ([Smith fc Terrild 
19841: iParesce fc Burrowsl 119871: iGolimowski et al l [ 1993 1 : 
Lecavelier des Etangs et all 119931: IKalas fc Jewittj | 1995[ 
Mouillet et al.lll997t iHeap et al.ll2000l: IKalas et allhoool ; 
Golimowski et al.l [2006); (2) imaging observations of 



dust th ermal emissio n in th e mid-IR (Lagag e~fc Pantinl 
| 1994l: iPantin et all 119971: iHeinrichsen et alj 119991: 
IWahhaj et all 120031: ITelesco et al.l l2005h and sub- 
millimeter ([Holland et all I1998D: (3 ) mid-R narrow 



band photometry (|Telesco fc Knackd I1991D. as well as 
broadband photome try in the mid- IR dGillettl 19861 : 
| Backman et~aIlll992D and millimet er (IChini et all 1 19911: 
Liseau et al.1 |2003[) : (4) optical dGledhill et all Il99ll: 
| Wolstencroft et all 119951 ) and near-IR (|Tamura et al.l 
l2006h imagin g of polarized ligh t ; and (5) mid-IR 



spectroscopy (Ajbken et al.l 119931: IKnacke et al.l 119931 : 



IWeinberger et al.ll2003l : lOkamoto et al. I I2004D . 

To accurately describe the dust grains in the (3 Pic- 
toris disk, all of the aforementioned observational data 
must be modeled simultaneously. While recognizing the 
importance of a comprehensive, simultaneous modeling 
of all of the j3 Pictoris observations, in this paper we 
present a simple model focusing on our Spitzer IRS spec- 
trum. We defer the comprehensive modeling of the dust 
around (3 Pictoris that attempts to simultaneously fit all 
of the observed data to a separate paper (A. Li et al. 
2007, in preparation). 

We model the Spitzer IRS spectrum of Pictoris 
building on the work of ILi fc Greenberg! (jl998l ) which 
includes (1) a population of cometary dust, composed 
of low-density porous aggregates of amorphous silicate 
core-carbonaceous mantle grains, 12 and (2) a population 
of crystalline olivine aggregates with the same size and 
spatial distributions as the amorphous cometary grains 
(but actually we do not require crystalline dust outside 
~ 60 AU from the star since at distances > 60 AU from 
the star silicate grains are too cold to emit at the char- 
acteristic mid-IR bands). 

For the amorphous cometary and crystalline olivine 

11 Indeed, as shown in Figure la, the SL point-source-extracted 
spectrum centered on the star (with a beam diameter of ~3.6", 
approximately the same as lKnacke et al. ( 1993)'s spectroscopic ob- 
servation) is close to that of [Knacke et al. (1993). 

12 The mass ratio of the silicate core to the carbon mantle is 
assumed to be m car t, /m s n 0.7 based on the cosmic abundance 
consideration (see Appendix A of Li & Lunine 2003). 
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fluffy aggregates, we assume a vacuum volume fraction 
P = 0.90 (see §2 of Li, Lunine, & Bendo 2003 for justifica- 
tion) and approximate the porous grains as spheres with 
radii a. We assume a power-law dust size distribution 
dn(a)/da oc a~", with a minimum grain radius a m j n 
1 /xm., and a maximum grain radius a max = 10000 /im 
(see Li, Luni ne, fc Bendo 2003 for ju s tificat ion). 

Similar to lArtvmowicz fc Clampinl (|1997l ). we assume 
a modified power-law spatial distribution for the dust, 

n(r) oc 1/ /(r/vp)^ 1 + (r/r p ) 7 . This functional form, 

peaking at r p x (//V) 7 , behaves like a power-law 
n(r) oc r -7 at large distances (r > r p ). At r < r p , 

the (r/r p ) 1 term dominates, approximating an increas- 
ing profile in the inner, dust-relatively-depleted region. 
We take r p = 100 A U and 7 = 2.7 as derived from scat- 
tered light modeling (|Artvmowicz et alll989f) . The inner 
boundary r m i n of the disk is taken to be the location in- 
side which silicate dust sublimates. For the (3 Pictoris 
disk, micron-sized silicate grains possess an equilibrium 
temperature T w 1500 K at r « 0.2 AU; therefore, we 
take r min = 0.2 AU. We take r max = 2000 AU. 

In modeling the dust IR emission, we - with all other 
parameters (i.e. r p , 7, r mm , r max ) pre-chosen- are there- 
fore left with only 2 free parameters: a - the dust size 
distribution power-index, and / which determines the 
amo unt of dust in the inner d i sk reg ion r < r p [we note 
that lArtvmowicz fc Clampinl (fl997) took /=1]. It is 
found that, with a w 3.2 and / ks 0.15, we are able 
to obtain a reasonably good fit to the overall SED and 
the Spitzer IRS spectrum including the silicate emission 
features at 10, 28 and 33.5 /im (see Figure |2J). 13 The 
required dust masses are approximately 1.63 x 10 27 g « 
0.27 M ffi and 4.89 x 10 25 g rj 0.0082 M e for the amor- 
phous and crystalline components, respectively. 

As shown in Figure [2}d, the model is somewhat de- 
ficient at A ~5-8/im and A ~ 17-24 /im. This is due 
to the oversimplified dust spatial distribution function 
adopted above. By including 2 rings or clumps of dust 
at r < 10 AU an d at r ~ 20-30 AU as im plied by 
iTelesco et ail (|2005f ) and lOkamoto et al.l (|2004[) . it is ex- 
pected that the model will account for the excess emis- 
sion at 5-8 /im and 17-24 /im (see A. Li et al. 2007, 
in preparation). The oversimplification of the assumed 
dust spatial distribution is also reflected by the model- 
predicted IR emission from the dust within a 3.6" di- 
ameter (corresponding to ~35AU; see Figure 2a): the 
model-predicted emission from the inner 35 AU region ac- 
counts for almost all the observed A < 15 /xm emission for 
the entire disk. Although this is consistent with the pre- 
vious spatially resolved mid-IR spectroscopy which show 
that 10 /Ltm silicate emission originates from the inner 

13 The amount of dust in the inner disk (r < r p ) is inverse 
proportional to /. With / = 0.15 and 7 = 2.7, the dust spatial 
distribution peaks at 0.46 r p 46 AU. At a first glance, this ap- 
pears inconsistent with the scattered light images which suggest 
a dust spatial distribution slo pe change at ~ 80-120 AU (e.g. see 
Golimowski et al. (1993, 2006)). However, the dust spatial distri- 
bution from 46 AU (~ ^p) to r p is rather flat: n(r) only changes 
by a factor of f»2.3 from ^r p to r p , while it drops much more 
steeply from r p to 2r p (by a factor of Ri5.8). In a subsequent pa- 
per (A. Li et al. 2007, in preparation), we will investigate in detail 
whether the present dust spatial distribution is able to reproduce 
the scattered light images. 



~20AU (jWeinberger et all 120031: lOkamoto et al.ll2004h . 

it appears to contradict the IRS detection of an apprecia- 
ble amount of silicate emission from the outer ~ 35 AU 
(see Figure la). This will be investigated in detail in a 
subsequent paper (see A. Li et al. 2007, in preparation). 

4. GAS MASS UPPER LIMITS 

Core-accretion models suggest that giant planets ac- 
crete their gaseous envelopes on timescales between 1 
and 20Myr. Therefore the (3 Pictoris disk, with an age 
of ^12Myr, is an excellent source to search for bulk 
gas. If the gas has a solar composition, then the bulk 
gas is expected to be hydrogen. ISO observations indi- 
cated line emission from H2 S(0) at 28.2 fim and S(l) 
at 17.0 /mi, suggesting that the disk possesses 54 M© 
warm H2 with an excitation temperature T ex = 110K 
(|Thi et al.ll200l"l ). However, apparently conflicting FUSE 
observations constrain the circumstellar H2 Lyman se- 
ries absorption and place a 3er upper limit of A^H^) < 
10 18 cm~ 2 on the H2 column density, significantly lower 
than (5-500) x 10 20 cm" 2 expected if the ISO-detected 
H2 we re uniformly distributed in an edg e-on disk in the 
beam (|Lecavelier des Etangs et al.l2001l ). an assumed ge- 
ometry that is consistent with observations of the atomic 
gas. Detailed studies of gas drag on the dust dynamics 
suggest that the gas:dust ratio is less than 1 (or <0.4 
molecular hydrogen exists in the disk). For exam- 
ple, if the disk possessed 40 M e gas, then small grains 
would collect at distances >200 AU and would increase 
the scatte red light surface brightness b y more than a fac- 
tor of 10 (jThebault fc Augereaull2005h . Searches for H I 
21 cm emission constrain the mass of atomic hydrogen 
Mhi < 0.5 M e (|Freudling et al.l [l995h : recent chemical 
models constrain the total mass of hydrogen <15 M ffl 
including t he molecular com ponent based on these ob- 
servations (|Kamp et al.ll2007| ). 

We searched for emission from H 2 S(2), S(l), S(0) and 
[S I], Fe II, and Si II but did not detect any of these 
species. We place 3<r upper limits on their line fluxes 
toward (3 Pictoris (see Table 1). Our 3<r upper limit on 
the H 2 S(l) line flux is <1.2xl0 _14 ergs _1 cm -2 , a factor 
of ~ 6.4 times lower than the reported ISO detection (see 
Figure [3J . Although the Spitzer IRS SH slit is a factor of 
two shorter than the ISO SWS slit (14"x27"), our upper 
limits effectively constrain the H2 line emission because 
the warm, bulk H2 is expected to be located at radii 
<100AU(<6"). 

Converting 3er upper limits on the line fluxes from any 
species into gas masses depends sensitively on the as- 
sumed gas temperature. Detailed models of the heating 
and cooling of molecular and atomic gas via gas-grain col- 
lisions, cosmic rays, line emission, etc. have been used 
to infer the temperature and chemical structure of the 
(3 Pictoris disk. However, these models depend on the 
initial gas:dust ratio or gas mass assumed. For exam- 
ple, a 2 M© disk (with an interstellar gas:dust ratio of 
~100) heated primarily by gas-grain collisions may pos- 
sess gas as warm as 100-150 K at intermediate heights at 
distances of 300-500 AU where [O I] fine structure emis- 
sion at 63.2 /im is the dominant coolant; while, a 0.2 M© 
disk may possess gas as warm as 300 K at lower heights 
and similar radii in the disk where [O I] and H2 rota- 
tional/vibrational emission a re the dominant coolants 
(|Kamp fc van Z adclhoff I2001D . Estimates of the total 
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gass mass, inferred from the measured column densi- 
ties of atomic species and the scattered light gas den- 
sity profile of Na I, suggest that the (5 Pictoris disk con- 
tains ~7.4xl0 -4 M® measured gas, corresponding to a 
gas:dust ratio ~0.019, significantly less than assumed in 
many detailed chemical models. Even if the gas has a so- 
lar hydrogen abundance relative to the heavy elements, 
then the disk possesses ~3.7xl0 -3 M® gas, correspond- 
ing to a gas:dust ratio ~0.093 (A. Roberge, private com- 
munication), still significantly less than the gas:dust ra- 
tios of 100 or 10 assumed in these models. 

To estimate an upper limit on the mass of H2 based on 
our S(l) line flux upper limit, we must assume a gas ex- 
citation temperature. Gas temperatures can be inferred 
from (1) detailed models of gas in thermal balance that 
calculate the composition and density structure of disks 
in addition to the temperature structure (as described 
above) or (2) observations if multiple electronic transi- 
tions are observed. Far-UV observations of f3 Pictoris 
have constrained the bulk circumstellar gas temperature. 
Analysis of the C I ( 3 P) A1613 and A1561 multiplets 
suggest that the excitation t emperature of the sta ble C I 
component, T cx = 50-100 K (jRoberge et afll200tf ). Since 
bulk hydrogen with a gas:dust ratio of 100 has not been 
directly detected (as assumed for the Kamp & van Zadel- 
hoff models) , we rely on the measured C I excitation tem- 
perature to act as a guide for the bulk gas kinetic tem- 
perature Tki n = 50-100 K. Therefore, we place 3<r upper 
limits on the mass in each of the listed species assuming 
gas temperatures T cx = 50 and 100 K (see Table 1). The 
total flux produced by N atoms or molecules 



F 



III 



(1) 



where d is the distance to the star, E = hv is the energy 
of the radiated photons, \ u is the fraction of atoms or 
molecules in level u, and A u i is the transition probability. 
Since the temperature of any H 2 is uncertain, we plot the 
3<7 upper limits on the H2 mass obtained from our S(l) 
and S(0) line fluxes as a function of temperature (see 
Figure 0J. If T gas = 100K, we estimate that <17M® H 2 
hydrogen remains in the disk from our H2 S(l) upper 
limit. 

The nondetection of H 2 , [S I], Fe II, and Si II emission 
toward /3 Pictoris is consistent with Spitzer observations 
of other 10-30 Myr old disks. ISO observations of the 
20Myr old dusty A1V star 49 Ceti have indicated line 
emission from H2 S(0) at 28.2 /im, suggesting that the 
49 Cet disk possesses >110M® warm H 2 with a n exci- 
tation temperature T ox < 100 K (|Thi et al.ll200lh . How- 
ever, high resolution Spitzer IRS observations place 3cr 
upper limits on the H 2 S(0) line flu x that are a facto r 
~9 lower than the reported detection (|Chen et al.ll2006l) . 
A recent Spitzer IRS search for H 2 emission from warm 
gas and millimeter search for 12 CO from cool gas around 
8 sun-like systems with ages <30 Myr measured upper 
limits <0.04 Mj up gas within a few AU of the disk in- 
ner radius (1-40 AU ) and <0.04 M® gas at 10-40 AU 
(jPascucci et al.l[2006h . Spitzer IRS observations place 5cr 
upper limits on the line flux from H 2 S(l) and S(0) to- 
ward HD 105, a GOV member of the Tucana-Horolgium 
association at 40 pc from the Sun, suggesting that less 
than 1500, 12, and 0.95 M®, respectively, at T gas = 50, 
100, or 200 K remains in this disk at an age of 30 Myr 



(|Hollenbach et al. 2005). Together, these observations 
suggest that warm molecular hydrogen, with T ex = 100 
K, dissipates from circumstellar disks on timescales <10 
Myr. 

5. THE ORIGIN OF THE ATOMIC GAS 

The origin of the observed stable atomic gas in Keple- 
rian rotation around (3 Pictoris (described in the Intro- 
duction), like the micron-sized dust grains, is currently 
not well-constrained. There have been suggestions in 
the literature that the stable g as component is produ ced 
by infalling refractory bodies dBeust fe Valironll2007l) or 
by collisions between dust grains (jCzechowski fe Mannl 
l2007t iFernandez et al.l l2006) . We consider these possibil- 
ities and hypothesize that some of the observed atomic 
sodium is generated by photon-stimulated desorption 
from Si0 2 -like surfaces. 

Infalling refractory bodies undoubtably contribute at 
least some portion of the observed gas but whether all of 
the gas is generated by the evaporation of parent bodies 
is uncertain. The spatial distribution of the gas may pro- 
vide a clue to its origin. The atomic gas is detected to 
distances of 150-200 AU on the southwest side and dis- 
tances of 300-350 AU on the northeast side of the disk. 
These distances correspond to the inner regions of the 
dusty disk that has been detected in thermal emission 
and scattered light to distances >1400 AU. The SW/NE 
dust brightness asymmetry observed in scattered light is 
also seen in the resonantly scattered gas emission lines; 
however, the asymmetry in the Na I emission is more 
pronounced. Finally, the gas and dust populations share 
an in ner disk that is tilted ~ 5° with respect to the outer 
disk (|Brandeker et al.| [2004). If infalling refractory bod- 
ies produce all of the observed gas, why are the spatial 
distributions of the gas and dust so similar? 

To explain the similarity in the s patial distributions of 
the du st and gas around (3 Pictoris. ICzechowski fe Marml 
(2007) recently proposed that the atomic gas is pro- 
duced via collisions between dust grains that vapor- 
ize at least some portion of the dust. Vaporization of 
dust grains in energetic collisions is believed to pro- 
duce atomic gas in supernovae shocks and may gener- 
ate a non-trivial fraction of the gas observed in the ex - 
osphe res of Mercury and the mo on dKillen et al.ll200lh . 
In the lCzechowski fe Mannl (|2007h model, radiation pres- 
sure accelerates sub-blowout sized grains ((3 meteoroids) 
to large radial velocities; these small particles then col- 
lide with bound, orbiting grains with relative veloci- 
ties as high as 90 km/sec (for carbonaceous grains). 
The observed abundances of the f3 Pictoris dust and 
atomic gas may provide a clue to the origin of the stable 
atomic gas. If the gas is liberated in vaporizing collisions 
and no atomic species are selectively removed (e.g. via 
out gassing or differentiation), then the gas composition 
should be the same as the grain composition. Our mod- 
els for the dust grains around f3 Pictoris reproduce the 
Spitzer IRS spectrum assuming a carbon relative to sil- 
icon abundance ratio of ^8.9, significantly smaller than 
the measured ratio of the c olumn densities of ato mic car- 
bon and silicon gas, ~500 (jRoberge et al.ll2006l ). 

Photon-stimulated desorption (PSD) may produce a 
substantial fraction of the sodium observ ed in the tenu- 
ous exospheres of Mercury and the Moon dLeblanc et al] 
120061 iLammer et all 120031: lYakshinskiv fe Madevi 11999?) 
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and may produce some of the spatially-extended sodium 
gas observed via resonant scattering toward Pic- 
toris. Since PSD produces gas directly from the 
dust, it might naturally explain the similarity in 
the gas and dust density distributions. The cross- 
section for PSD jumps dramatically at UV wave- 
lengths, becoming most e fficient at A < 2500 A 
(|Yakshinskiv fc Madevi fl999l ) . Pictoris has been ob- 
served extensively in the far and near ul t ra violet using 
FUSE (|Lecavelier des Etangs et alj|2001t iRqberge et al'l 
[20061 ). HST GHRS and STIS ( |Roberge et al.l l2000l). with 
a continuum flux at 2000 A < A < 2500 A measured 

; A _1 G 



by I U E, F x > 4 x IP" 11 erg s" 1 cm 



janz 



et al.l 

Il995h . iFernandez et all (|2006l ) model the observed stellar 
ultra violet spectrum of Pic using a rotationally broad- 
ened (v Iot = 130±4kmsec" 1 ) PHOENIX model with ef- 
fective temperature T* = 8000 if, and surface gravity, 
logg = 4.2, suggesting a stellar UV photon production 
rate, L > = 2.0xl0 44 s _1 . 

' 1000-2500 A 

The efficiency with which photons desorb sodium, av- 
eraged over the surface of a sphere, is 



(2) 



(jYakshinskiv fc Madevi 1 1999T ) where Q is the photon- 
stimulated desorption cross section and a is the atomic 
surface coverage. Laboratory experiments measure Q^a 
~ (3±1) x 10~ 20 cm 2 for lunar samples with tempera- 
tures of 250 K and using incident photons with A < 2500 
A, suggesting an efficiency of £Na,250K = 2.3 x 10~ 8 for 
materials with a lunar temperatures and compos itions, 
a Na = 3 x 10 12 cm" 2 (jYakshinskiv fc Madevi fl999[ ). In 
addition, the photon-stimulated desoprtion efficiency of 
sodium is also function of substrate temperature, T gr . 
We estimate 



e Na = (6.93 x 10~ 8 )10- 



-122.8K/T gr 



(3) 



from measu rements of sodium desorpti on yields in the 
laboratory fYakshinskiv & Madevi [2004D . normalized to 
the sodium desorption efficiency at the moon as discussed 
above. 

The rate at which atomic sodium is produced in the 
Pictoris disk is 



lixrdr 



L 



where r = 120 AU, W = 6.6 (r/r ) a75 , and r D = 2x 10~ 3 . 
We assume that the bulk of the surface area is con- 
tained in 1 /im amorphous olivine grains and estimate 
the grain temperature as a function of distance using Mie 
theory (the grains are spherical ) and using laboratory - 
measured indices of refraction (jDorschner et al.l fl995). 
If R* = 1.69 R and T* = 8000 K, we estimate a 
sodium production rate 1.3 x 10 33 s _1 from the grains 
via photon-stimulated desorption. We do not use the 
same dust spatial distribution as adopted in §3 for mod- 
eling the IR emission because the desorption properties 
of porous aggregates invoked in §3 are poorly known, 
while the desorption properties of compact grains (as- 
sumed in eq.[4])are better quantified (e.g. experimental 
data are available). Moreover, the dust spatial distribu- 
tion adopted in §3 is already oversimplified. 

We estimate the number of sodium atoms in the Pic- 
toris disk from observations of resonantly scattered Na I 
and infer the time required to produce the circumstellar 
sodium. Since the ionization potential of sodium is 5.1 
eV and Pic possesses a high UV luminosity, the major- 
ity of the circum s tellar sodium is expected to be ionized. 
Fernandez et al.l (|2006[ ) have written a photoionization 
code that calculates the densities of neutral and ionized 
elements from H to Ni, assuming that the Pic disk is 
optically thin and has a solar composition; they estimate 
a sodium neutral fraction of 3xl0 -4 at the disk mid- 
plane at 100 AU from the central star. Since the sodium 
neutral fraction is a function of distance from the stellar 
ionizing source, we have written a simple photoionization 
code which calculates the densities of neutral and ionized 
carbon and sodium in the disk as a function of position 
assuming that (1) the gas is optically thin, (2) the disk 
gas contains only carbon and sodium, (3) ionization of 
carbon produces the electron density in the disk, (4) the 
gas temperature, T cx — T gr , and (5) the circumstellar 
Na I possesses a density distribution 



n-Nai = n a 



-1/2 



exp 



(-)' 



(6) 



dz ■ 



1000-2500 A 



epdust(r)] rra 2 n dus t 



47rr 2 

(4) 

where i? m i n = 3AU and -R max = 232 AU are the min- 
imum and maximum radii, and z m j n = — 100AU and 
z max = 100 AU are the minimum and maximum heights 
at which the atomic sodium gas is observed. Since the 
Na I gas in the Pic disk is located at distances of up 
350 AU, well beyond our IRS slit, we use a dust density 
distribution inferred based on observations of the disk on 
large scales (rather than the dust density profile that we 
derive in Section 3). Artymowicz (private communica- 
tion) has inferred the quantity 7ra 2 ndust, s ummed over 
all gr ain sizes, from HST STIS observations (|Heap et all 
2000) assuming a grain albedo 0.5: 



Tra n dust 



W 



-1/2 



exp 



(w) 



0.7 



(5) 



([Brandeker et alj [2004) where n = (1.02 ± 0.04) x 
10- 5 cm- 3 , r = 117 ± 3AU, a = 0.94 ± 0.06, b = 
6.32 ± 0.04, and a = 0.168 ± 0.05. We use pub- 
^ Hsned carbon and s odium photoionization cross sections 
(IVerner et al.|[l996|) and radiative recombin ation coeffi- 
cients (jVerner fc Ferlandll 19961: iNahar | [l995) to estimate 
the carbon and sodium densities assuming that the atoms 
are in ionization equilibrium at each point in the disk. 
If the gas has a solar carbomsodium abundance ratio, 
then the Pictoris disk possesses 4.6 xlO 44 sodium atoms 
and a time of ~13000 years is required to produce the 
sodium gas. Models of the breaking of atomic gas in the 
Pic disk suggest tha t the lifetime of an atom in the disk 
is ~10 4 to 10 5 years ([Fernandez et al.|[2006h . suggesting 
that the sodium gas is in a steady state if it is produced 
by PSD and braked by Coulomb interactions. 

We can estimate the radial dependence of the sodium 
number density and compare it to that i nferred from res- 
onan tly scattered sodium observations (fBrandckc r et all 
I2004D . If atomic sodium produced in the disk does not 
migrate radially, then sodium gas generated by photon- 
stimulated desorption will have a number density distri- 
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bution 

L 

( \ 1000-2500 A \rr\ ( \\ 2 / \ 

nir, z) oc — — e Na [T dnat {r, z)\ ttcl n dust (r, z) 

(7) 

where r is the distance to the central star and irc^ridust 
is the surface area contained in the dust grains. The 
predicted PSD-produced sodium density distribution at 
the disk midplane falls too quickly as a function of ra- 
dius compared to that inferred from observations (see 
Figure 5). The rapid fall off in predicted sodium density 
is the result of the 1/r 2 dilution of stellar photons. The 
sodium desorption efficiency is able to partially overcome 
this effect because eNa[Tdust{r, z)] oc r if the grains are 

1 /im olivine spheres. However, even taking into account 
the temperature dependence of the desorption efficiency, 
our predicted sodium number density falls off too quickly 
with radius by a factor of r. If PSD produces the bulk 
of the observed sodium, sodium atoms must migrate to 
larger distances via radiation pressure. 

Our model for the origin of circumstellar sodium could 
be extended to other atomic species in the (3 Pictoris 
disk. At the current time, few laboratory measurements 
have been made for the photon-stimulated desorption 
rates of other atoms. Potassium is the only species, 
other t han sodium, for whic h PSD rates have been mea- 
sured. iMadev et al.l (jT998) measure cross sections for 
postassium desorption from Cr203 surfaces, Qk — 2 x 
1CT 20 cm 2 at 2500 - 3500 A and 2 x 10" 19 cm 2 at 1900 
A. If photon-stimulated desorption of potassium has the 
same temperature dependence as observed for sodium 
and the (3 Pic circumstellar grains have a solar potassium 
to sodium abundance ratio, ok = 2.1 x 10 11 cm~ 2 , then 
we estimate a potassium production rate in the (3 Pic 
disk of 3.3 x 10 32 K atoms s -1 , ~3000 times less than 
the sodium production rate. Atomic potassium has not 
yet been detected in the /3 Pictoris disk and 3(7 upper 
limits on the circumstellar potassium mass have not yet 
been reported in the literature. 

6. DISCUSSION 

The simple dust model presented in §3, consisting of 
cometary grains with a power-law size distribution and a 
modified power-law spatial distribution, provides a rea- 
sonably good fit to the overall SED and the observed 
Spitzer IRS spectrum. At a first glance, the model has 
many free parameters: (1) a m ; n , a max and a for the dust 
size distribution dn/da; and (2) r m - m , r max , /, r p , 7 for 
the dust spatial distribution dn/dr. With general con- 
straints from the disk structure and the dust absorption 
and emission properties, we actually have only 2 free pa- 
rameters a and / left (see §3 for details). In so doing, 

2 assumptions have been made: (1) all grains have the 
same spatial distribution; and (2) the grains producing 
the IRS emission spectrum are also responsible for the 
optical scattered light or their spatial distribution fol- 
lows that of the dust responsible for the optical scattered 
light. Admittedly, the former assumption is oversimpli- 
fied since the response of the grains to the stellar radi- 
ation pressure and Poynting-Robertson drag varies with 
grain size and therefore the simple radial power law is 
probably not all that representative of the actual grain 
distribution uniformly with size. To justify the latter as- 
sumption, a simultaneous modeling of both the scattered 



light and IR emission is required. Nevertheless, the gen- 
eral conclusion for the dust (the presence of crystalline 
silicates and an inner warm dust component in the disk) 
remains valid. 

We have presented a model for the production of cir- 
cumstellar sodium gas around (3 Pictoris from dust via 
photon-stimulated desorption. Our model is able to gen- 
erate the inferred circumstellar sodium within the res- 
idence time of the gas but is unable to reproduce the 
radial dependence of the number density distribution if 
the sodium atoms do not migrate to larger distances via 
radiation pressure. Our estimate for the time required to 
desorb the observed sodium is a function of (1) the cir- 
cumstellar grain temperature, (2) the circumstellar gas 
temperature, and (3) the gas carbomsodium abundance 
ratio. 

The photon-stimulated desorption rate of sodium from 
orbiting dust grains is dependent on the grain tempera- 
ture with sodium more efficiently desorbed from warmer 
surfaces. Since the bulk of the surface area of the dust is 
contained in the smallest grains, we estimated the grain 
temperature assuming that the grains have radii of a = 1 
/im, are composed of amorphous silicates, and are spher- 
ical. However, if the grains are very small (27ra < A) 

then the grain temperature, T gr = [0.25(i?*/r) 2 ] °' 2 T*, 
and the sodium desorption rate, T = 3.7 x 10 33 s" 1 , and 
4000 years are required to produce the sodium gas. Or, 
if the grains are very large (black bodies), then the grain 
temperature, T gr = [0.25(i?*/r)] ' T*, and the sodium 
desorption rate, T = 2.5 x 10 32 s _1 , and 58000 years are 
required to produce the sodium gas. 

In our simple model, we have assumed that the gas 
and dust have the same temperatures and we have ex- 
plored the possible range in dust temperature profiles. 
We also consider how changes in the gas temperature 
profile affect our estimates for the sodium neutral frac- 
tion and therefore the total number of sodium atoms in 
the disk. The estimated neutral fraction of sodium is 
weakly dependent on the gas temperature and composi- 
tion. The radiative recombination coefficients for carbon 
and sodium decrease by factors of ~5 if the gas temper- 
ature increases from 10 K to 100 K or 100 K to 1000 
K. If the gas has a temperature profile similar to small 
grains or black bodies (instead of 1 /im olivine spheres), 
then we estimate that the (3 Pictoris disk possesses cither 
7.0xl0 44 or 5.4xl0 44 sodium atoms, respectively, similar 
to the 4.6 x 10 44 sodium atoms expected if the gas has the 
temperature profile expected from 1 /im amorphous sili- 
cate grains. The assumed abundance of sodium relative 
to carbon also affects the total number of sodium atoms 
inferred from the neutral sodium observations. In our 
model, we have assumed that the gas has a solar compo- 
sition; however, a recent inventory of the circumstellar 
atomic gas around (3 Pictoris suggests that carbon is en- 
riched compared to iron and oxygen b y factors of 16 and 
18, respectively (jRoberge et al.l [20061 ). If the j3 Pictoris 
disk is enriched in carbon relative to sodium by a factor 
of 16, then the disk possesses 1.7xl0 44 sodium atoms, 
somewhat less than the number inferred if the gas has a 
solar composition. 



7. CONCLUSIONS 
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We have obtained Spitzer Space Telescope IRS spectra 
of the 12 Myr old debris disk around Pictoris. We find 
that: 

1. In addition to the 10 /xm silicate emission feature 
originally detected in ground-based observations, we, for 
the first time, also observe weak crystalline silicate emis- 
sion features at 28 /zm and 33.5 /xm. 

2. The IRS dust emission spectrum and the overall 
SED are well reproduced by a dust model consisting of 
fluffy cometary and crystalline olivine aggregates. 

3. No H 2 , [S I], Fe II, or Si II emission is detected. Our 
3cr upper limits suggest that <17M® H2 remains in the 
disk significantly less than the previously reported de- 
tections of H2 S(l) and S(0) emission. The circumstellar 
disk around Pictoris has too little gas to support the 
formation of giant planets. 

4. Some of the observed resonantly scattered Na I, 
observed at visual wavelengths, may be produced via 
photon-stimulated desorption; the timescale to generate 
all of the inferred circumstellar sodium is 13000 yr, ap- 
proximately the residence time of the gas in the disk if 
the gas is braked by Coulomb forces. 
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TABLE 1 

Spitzer IRS Line Flux Upper Limits for (3 Pictoris 









50 K 


100 K 


Atomic/Molecular 


Gas Species 


Wavelength 


Line Flux 


M gas 


Mgas 


Data 




Om) 


( ergs -1 cm -2 ) 


(M e ) 


(M e ) 


References 


H 2 S(2) 


12.279 


<5.0xl0~ 14 


<1.6xlO n 


<1.7xl0 4 


2,4 


H 2 S(l) 


17.035 


<1.2xl0- 14 


<2.1xl0 5 


<17 


2,4 


[S I] 


25.249 


<8.6xl0- 14 


<0.047 


<0.00016 


1 


Fell 


25.988 


<6.1xl0~ 14 


<0.29 


<0.0011 


3 


H 2 S(0) 


28.221 


<4.1xl0~ 14 


<3400 


<43 


2,4 


Si II 


34.814 


<4.9xl0~ 13 


<0.028 


<0.00046 


1 



References. — (1) Haas, Hollenbach, & Erickson 1991; (2) Jennings, Weber & Brault 1987; (3) Nussbaumer & Storey 1980; 
(4) Wolniewicz, Simbotin, & Dalgarno 1998 
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Fig. 1. — (a) Spitzer IRS SL, SH, and LH spectrum of j3 Pictoris (red line) with the PHOENIX atmospheric spectrum (black line) 
overlaid. The inserted panel illustrates the 5-36 fim stellar-subtracted Spitzer IRS du s t exce ss emission. Also plotted are the UBVRIJHK 
stellar photometry (magenta squares), the ground-based spectrometry of[Knackc ct al. ( 1993) (lab eled " K93"), and the IRS spectrum in the 
inner 3.6" [corresponding to ~35AU; approximately the same beam size as that ofKnackc ct al. (1993)] disk (green line), (b) Comparison 
of the stellar-subtracted Spitzer IRS dust excess emission (red line) with the IRAS photometry (cyan open circles), the Big Mac/IRTF 
p hotometry (black fille d circles), the IRTF N band and Q band photometry (green open stars), and the IRTF 2.6-13.5 fim spectrometry 
of lKnacke et al.l | |1993D (blue open squares). Numerous hydrogen absorption features can be seen between 5 and 13 /im; their depths are 
consistent with that expected for a star with A5V photosphere at spectral resolution of the IRS. 
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Fig. 2. — Comparison of the model (black line) to (a) the observed dust IR emission spectral energy distribution and (b) the Spitzer 
IRS spectrum ( red line). Also plott ed is the model-predicted emission in the inner 35 AU disk. Blue open squares: the IRTF 2.6-13.5 ^m 
spectrometry of Knackc ct al. (1993) increased by a factor of 1.94. Black filled circles: the Big Mac/IRTF photometry of[T clcsco & Knackc 
H1991I1 increased by a factor of 1.94. Green open stars: the IRTF N band and Q b and ph o tomet ry of IBackman et al.l (|1992l) with the 
N band incre ased by a factor of 1 .94. Cyan open circles: the IRAS photometry of Gillctt (1986). Magenta filled triangles: the KAO 
photometry of[Harvcv ct al. (1998). Black filled squares: the 850 /im SCUBA photometry of Holland ct al. (1998). With this increase (of 
a factor of 1.94), all observational data agree with each other very well. The inserted panel in (b) illustrates the model fit to the 10 /im 
silicate features. 
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Fig. 3. — (a) Portion of the observed Spitzer SH spectrum of j3 Pictoris around H2 S(l) with a model for unresolved H2 emission 
(overplotted in red) inferred from ISO detection, (b) Same as (a) but for the Spitzer LH spectrum around H2 S(0). 
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Fig. 4. — Three a upper limits on the hydrogen gas mass as a function of disk gas temperature, estimated from the line fluxes of H2 
S(l) (solid line), and H 2 S(0) (dotted line). 
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Fig. 5. — The estimate d number density dist ribution of Na I atoms at the midplane of the j3 Pic disk, inferred from observations of 
resonantly scattered Na I (Brandekcr ct al. 2004), is shown with triangles. For comparison, we overplot the number density distributions 
of Na gas produced via PSD, assuming that the grains are small (2na < A; dotted line), large (black bodies; dashed line), or 1 fim olivine 
spheres (dashed-dotted line). Our models are normalized such that they contain the same mass as the fitted model. 



